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I. INTRODUCTION 



Recently, there has been a renewed interest in the investigation of the dispersion forces 
which originate from the existence of zero-point electromagnetic oscillations. The most well 
knownphysical phenomena of this kind are the van der Waals and Casimir forces (see, e.g., 



Refs 



.EH. 



Both forces are of the same nature and are related to nonrelativistic and 



relativistic situation, respectively. In the last few 



force between metals were performed 



years many measurements of the Casimir 



Q, B, 0, B, Q, E3, EH E3, Q Q 



The increasing interest 



in experimental research of dispersion forces is motivated by their role in nanotechnology 



16 
14 



11 

17j and for constraining predictions of modern theories of fundamental interactions lc 

~nn n 

18U19U20H (see also Ref. |2jj covering both the experimental and theoretical developments 



on the subject ). Quit e recently both the van der Waals and Casimir forces were shown to be 
imp o rtMt fiHQ in expend on o,an tum refl ec tl0 n aa d Bose-Eiastein condensation. 

The theoretical description of the Casimir force between real metals at nonzero temper- 
atnre is prohle-natic. The d lr ee t apphcation Q H Q of the Lhshi, tannla Q 
in combination with the dielectric function of the Drude model leads 301 to a violation of 



Nernst's heat theorem in the case of perfect crystals, i.e., to a nonzero value of the Casimir 
entropy at zero temperature depending on the separation distance between the plates. As 
was shown in Refs. 3>ElI|> for crystals with defects or impurities Nernst's heat theorem is 
satisfied, so that the entropy of a fluctuating field is equal to zero at zero temperature. This, 
however, does not solve the problem since a perfect crystal is a truly equilibrium system with 
non-degenerate dynamic state of lowest energy. According to quantum statistical physics, 
the Nernst heat theorem is valid for any system possessing this property. It follows that any 
formalism applied to a perfect crystal must satisfy the Nernst heat theorem, whereas the 



approach of Refs. [2 



26 



27 



281 ] does not. It is notable that this approach predicts large 



thermal corrections to the Casimir force at short separations as compared with the case of 
ideal metals. 

Other theoretical descriptions of the thermal Casimir force between real metals are based 
on the use of the Lifshitz formula combined with the plasma model dielectric function 



32 



33 



34j or with the Leontovich surface impedance 30|, |35|. In the framework of these 



approaches the Nernst heat theorem is satisfied for both perfect lattice as well as lattice 
with some small concentration of impurities. They predict small thermal corrections to the 



2 



Casimir force at short separations in qualitative agreement with the case of ideal metals. 

Recent experiments reveal the possibility to test the validity of different theoretical ap- 
proaches to the thermal Casimir force. Thus, the first modern measurement of the Casimir 



force by means o 
Refs. 



26 



27 



a torsion pendulum 5] was 
2SH and consistent with Refs. 




to be in disagreement with 
. The most precise and 



accurate modern experiment by means o 
eludes the theoretical approach of Refs. 
with the theoretical approaches of Refs. [30J, |32j, |33 



a micromechanical torsional oscillator 




13 



cx- 



at 99% confidence and is consistent 
3511 . Notice that these conclusions 



concerning the comparison between theory and experiment are not universally accepted (see 
the discussion on this subject in Refs. [38, 39]). 



27, ]2M is the non- 



A distinguishing feature of the theoretical approach of Refs. |2i 
monotonous dependence of the Casimir pressure and the free energy of a fluctuating field 
on the temperature. There are some temperature intervals where the magnitudes of the 
Casimir pressure and the free energy decrease with the increase of temperature. As a 
result, there are temperature intervals where the entropy of the fluctuating field takes neg- 
ative values. The same is valid for hypothetical (nonexistent) dielectrics with a frequency- 
independent dielectric permittivity e = 100 or higher [27]. By contrast, in the approaches of 



Refs. 



3(1 32, 3 



34. 



35| for metals and also for some real dielectrics considered in Ref. 



the magnitudes of the Casimir free energy and pressure increase with temperature. There 
is, however, no general proof that the magnitudes of the Casimir free energy and pressure 
must increase with temperature for all real materials. 



As can be seen from Ref. 



221 ]. the magnitudes of the Casimir- Polder free energy and force 



for the case of an atom interacting with a metal wall may decrease within some temperature 
region. The respective entropy of the fluctuating field becomes negative. When taken into 



22 



from 



account that the Lifshitz formula for an atom near a metal wall is obtained 
the case of a rarefied dielectric wall spaced parallel to a metal wall, it may be supposed that 
the primary reasons for the nonmonotonous behavior of the free energy are inherent in the 
Casimir interaction between metal and dielectric. 

In this paper we investigate the Casimir pressure, free energy and entropy in a configu- 
ration of two thick parallel plates, one being metallic and the other dielectric. Previously, 
only two metallic or two dielectric plates were considered in the retarded regime (see also 



Refs. 



41 



42| where a so-called "unusual pair of plates" was treated - one of which is made 



3 



of an ideal metal and the other is infinitely permeable). In the nonretarded limit of the van 
der Waals force the Hamaker constants for the configuration of a metallic and a dielectric 
p,ate were canted in Re, Q. On the b a sis of the Lifshn, fonnu.a we nnd the ana.ytic 
expressions for the Casimir pressure, free energy and entropy in the configuration of a plate 
made of ideal metal and a parallel plate made of dilute dielectric separated by a gap of thick- 
ness a. It is shown that there are intervals on the temperature axis where the magnitude 
of the Casimir pressure decreases with the increase of temperature, whereas the magnitude 
of the free energy increases with the temperature and the Casimir entropy is positive. The 
asymptotic expressions of the obtained exact formulas at both low and high temperatures 
are presented and the validity of the Nernst heat theorem is demonstrated. The possibil- 
ity of a decreasing Casimir energy and of negative values of the entropy of a fluctuating 
field within some temperature intervals is demonstrated for the cases of ideal and real (Au) 
metals near dielectrics with arbitrary constant e (not dilute) and real dielectrics with the 
frequency-dependent dielectric permittivities (Si, a — Al 2 03). The physical interpretation 
of that negativeness of entropy is discussed. From this we give a positive answer to the 
above question of whether the entropy of a fluctuating field in the Casimir regime can be 
negative. This is important not only for the theory of the thermal Casimir force but also 
: or experimental investigations of the Casimir interaction between metal and semiconductor 
441- 

The paper is organized as follows. In Sec. II the main notation is introduced and the 
Casimir free energy, entropy and pressure are found analytically in the framework of the 
exactly solvable model (one plate is made of ideal metal and the other of dilute dielectric). 
Sec. Ill contains the results of analytical and numerical computations performed for the 
case of ideal-metal plate spaced parallel to the plate of dielectric with some frequency- 
independent e (not dilute). The case of Au plate and real dielectrics with the frequency- 
dependent dielectric permittivities (Si, a — AI2O3) is considered in Sec. IV. Sec. V contains 
our conclusions and discussion. 
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II. CASIMIR INTERACTION BETWEEN AN IDEAL METAL AND A DILUTE 
DIELECTRIC 



We consider the configuration of two thick parallel plates, one made of a real metal and 
another of dielectric at a separation a and temperature T in thermal equilibrium. The 



Casimir free energy per unit area is given by the Lifshitz formula 
dimensionless variables it can be represented in the form 



401 ] . In terms of 



^) = ^£(i-^.) 



1=0 

poo 

X 



ydy {In [l - rf (0, y)rjf(Ci, y)e~ v ] 
+ ln[l-r3 4 (C I ,yW(G,y)e-»]}. (1) 

Here, the dimensionless Matsubara frequencies are d = &/£ c , where £j = 2TtksTl/h are 
the usual Matsubara frequencies, and £ c = c/(2a) is the characteristic frequency, &b is 
the Boltzmann constant, 5/o is the Kronecker symbol. The reflection coefficients for two 
different polarizations of the electromagnetic field are expressed in terms of the dielectric 
permittivities e M,D (u) of a metal and a dielectric, respectively, 



M,D / 9 . /-it M,D 



ri (0,2/) = jl f = — , (2) 



where ' D = £ M ' D 

For later use we introduce the so-called effective temperature, T e g, defined through 
ksT e s = h£ c = he/ (2a). In terms of the effective temperature the nondimensional Mat- 
subara frequencies are expressed as d = tI, where r is a new parameter r = 2nT /T e Q. 

For an ideal metal rY ± {(i,y) = 1 and Eq. (JTJ) gives us the free energy of the Casimir 
interaction between an ideal metal and dielectric 



^) = ^£(i-^o) 

7— n V / 



1=0 

pea 

X 



ydy {In [l - rf(0, y)e^] 
+ ]n[l-r£(Ci,y)e-*]}. (3) 
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Let us assume that the dielectric is dilute, i.e., having dielectric permittivity ef = 1 + 77, 
where 77 = const <C 1. Expanding the reflection coefficients r® ± from Eq. (J2J) in powers of 
77, one obtains 



In [l-rf(Ci,y)e-y] = 77- 



4 V y 



(4) 



_ [(4 + e -v) r 4;4 _ 4 e -VZ 2 2/ 2 + 4 (e^ - 2) y 4 ] + 0(r/ 3 ), 



322/ 

In [l-r°(0,y)e- y ] = - V 



e-y r 2 l 2 

~1F 



Substituting Eq. (j3J) in Eq. Q and preserving only two powers in the small parameter 77, 
we arrive at the following expression for the Casimir free energy 



F m (a,T) 



her 

32vr 2 a 3 



E fi - 5*- 



Z=0 



x 



[4y 4 -e^ (r 4 / 4 -2r 2 /V + 2y 4 )] 



(5) 



Upon integrating and summing in Eq. (J3J), we obtain 



.P D (a,T) 



327r 2 a 3 



77T- 



e lT + 2re T - 1 
A(e T - l) 2 



(6) 



-77 



16 



2 (e 2r + 2re T 



e 4r + 4re 2T - 1 
4(e 



2r 



2 2r l -e 4T + 2r(l + 4e 2r + e 4r ) 

+T 6 2(e 2 - - l) 4 

00 00 
+2r 4 / 4 Ei(-2r/) - 2r 2 ^ Z 2 Ei(-2r/) 

1=1 1=1 

where Ei(z) is the integral exponent. 

The entropy of the Casimir interaction per unit area of the ideal metal and dilute dielectric 
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is found as minus the first derivative of Eq. © with respect to the temperature 



S LD (a,T) = 
A + t 



<9.P D (a,T) k 



dT 

1 ~r)e 7 



b_ l V 
8vra 2 | 2 



e T + 1 



-re 



-4re T . 



1+T 



1)3 

;i - r)e 



16 



2 r 



7e T + 9 



. 4 (e T - 1) 
2re 2T + e 2r - 1 



(7) 



(e T - l) 3 2(e 2T - l) 2 

3(1 - e 4T ) + 10r(l + 4e 2r + e 4r ) 



2(e 2 - - l) 4 

CO 

(6r 2 / 2 - 10r 4 / 4 ) Ei(-2r0 



i=l 



Eq. (jHJ) permits to find an analytic expression for the Casimir pressure between the ideal 
metal and dilute dielectric 

<9.F ID (a,T) Hct 



P ID (a,T) = - 

rV(e r + 1) 
+ (e--l) 3 



9a 
rf 

T 

4t 



32vr 2 a 4 
e 2T + 2re T - 1 



e 2r + 2re r - 1 
(e- - l) 2 



(8) 



r 2 e r (e T + 1) e 4r + 4re ir - 1 



8(e 



2t 



2(e 2r - l) 4 J J 



Note that Eqs. ©-© are perturbative in the small parameter 77 but for each perturbation 
order the thermal effects are taken into account exactly. This gives the possibility to find 
the high-temperature behavior (r > 1) of the Casimir free energy, entropy and pressure. 
Thus, from Eq. (jUJ) in the limit r ^> 1 (T ^> T e g) for the free energy it follows 



(9) 



In a similar manner from Eqs. (J7J) and (JBJ) one finds the behavior of the Casimir entropy and 
pressure at high temperature 



S m (a,T) 



did/ T ^ k sT ( 7 
P a , T) = -— — -7i 1 - —77 
lbna 6 \ 16 



(10) 



Needless to say that entropy and pressure from Eq. (|10|) can be obtained also as minus the 
derivatives of Eq. @ with respect to temperature and separation distance, respectively. 



Now let us consider the low temperature behavior of the Casimir free energy, entropy 
and pressure. For this purpose the direct application of Eqs. ©-© is rather cumbersome. 
The same results can be obtained more simply by using the Abel-Plana formula (see, e.g., 
Refs. flEil) 



E f 1 - 5*. 



1=0 



F(l) 



F(t)dt 



(11) 



+ i 



dt 



[F(it) - F(-it)} , 



where F(z) is an analytic function on the right half-plane. Before the application of Eq. ([lip, 
we return to Eq. (J3j) and perform the integration with respect to y (but not the summation 
over I). The result is 

her 



32vr 2 a 3 



E ( 1 - 5*. 



where 



F(l) 



[l + rl)e 



-rl 



V 
16 



4(l + rl)e- rl 



(12) 



(13) 



A (1 + 2rl + r 2 / 2 - 2r 3 / 3 ) e~ 2rl - 2r 2 / 2 (l - r 2 / 2 ) Ei(-2rZ) 

The sum in Eq. (|12|) can be calculated using the Abel-Plana formula (jllj) . In doing so, we 
note that the first term on the right-hand side of Eq. (jll|) . where F(t) is given by Eq. (|13|) 
with a substitution of t for Z, is proportional to the Casimir energy at zero temperature, 
whereas the second gives the thermal correction. Direct integration results in 



F(t)dt 



V 

T 



, 457 

1 77 

960 ' 



(14) 



The second integral on the right-hand side of Eq. (|TT|) is calculated as a perturbative ex- 
pansion in powers of a small parameter r with the aim to obtain the low temperature 
asymptotics. Preserving all terms up to r 3 inclusive, we obtain from Eq. (|13p 



F{it) - F(-it) = [r 3 t 3 + 0(r 5 )] 
3 

+ [_ 7rr 2 f 2 + 6r 3 f 3 + 0( ^ 



(15) 



Using this equation, the integral can be calculated 

dt 



TjT 



[F{it) - F(-it)] 



(16) 



?7 2 t 2 



720 32 



T 

To 



C(3) 



7T^ 



S 



where ((z) is the Riemann zeta function. 

Substituting Eqs. (fTT|). (JUJ), (JTfiJl in Eq. (|T2^) . the low temperature behavior of the free 
energy is obtained 

He 



f 1D (a,T) = - 
"457 



T 



32 



327T 2 a 3r/ i 14 720 
C(3)r 3 



(17) 



30 



7T Z 



T 

10 



[here the terms of order 0(t]t 6 ) and 0(v 2 t 5 ) are omitted]. 

The low temperature asymptotics of the Casimir entropy and pressure can be obtained 
by the calculation of minus the derivatives of Eq. ()12j) (with respect to temperature and 
separation, respectively) and by subsequent application of the Abel-Plana formula with 
the resulting functions F(l). The same asymptotic expressions are obtained also by direct 
differentiation of Eq. ([17)1 . They follows: 

k B 



S 1D {a,T) 
P ID (a,T) = - 



32na 2 ' 
he 



(jr rj 


[3C(3) t 


} 


\45 4 


_ 2tt 2 5_ 





(18) 



32vr 2 a 4 



T 

720 



320 



(457 



r 4 ) 



It is significant that S lD — > when T — > (recall that r ~ T by definition), so that the 
Nernst heat theorem is satisfied. For dilute dielectrics from Eqs. (J2J), p8)l it follows also 
that S lD > 0. At the same time, there are temperature intervals where the Casimir pressure 
is nonmonotonous. To demonstrate this, in Fig. 1 the relative thermal correction to the 
Casimir pressure, defined as 



5 T P 



P(a, T) — P(a, 0) 



(19) 



P(a,0) 

is plotted as a function of temperature for rj = 0.001 (solid line) and rj = 0.1 (dashed line) 
at a separation distance a = 2 /im [here we use Eq. (JHJ) and P(a, T) = P ID (a, T)]. As is seen 
from Fig. 1, the region where the thermal correction is negative is in accordance with the 



22| for the Casimir-Polder energy and force between an atom and a metal 



results of Ref. 
wall. 

From Fig. 1 it is seen that for rj — 0.1 it holds StP id < at temperatures T < 343 K. 
In this case he relative thermal correction has the minimum value StP id ~ —0.007 at 
T = 270 K. For comparison, at T = 400 K the thermal correction has the positive value 



9 



5 T P lD ~ 0.018. Note that there is some analogy between the behavior of the Casimir 
pressure and thermal correction as a function of temperature, obtained here on the one 



hand, and discussed in Refs. 

q y, 0, q 

on the other hand in the case of real metals 
described by the Drude dielectric function. We will return to this analogy in more detail in 
Sec. V after the extension of the above results for more realistic plate materials which will 
be obtained in the following sections. 

III. THERMAL CORRECTIONS TO THE CASIMIR INTERACTION BETWEEN 
AN IDEAL METAL AND DIELECTRIC WITH FREQUENCY-INDEPENDENT 
DIELECTRIC PERMITTIVITY 

In this section we consider one plate made of dielectric with arbitrary constant e D (not 
dilute) and the other plate made of ideal metal as before. The application regions of the 
asymptotic expressions at low and high temperatures, obtained in Sec. II, will be our initial 
concern. To gain an impression of how accurate the asymptotic of Eq. (fTTj) is, let us calculate 
the relative thermal correction to the Casimir energy, 

F(a,T) - F(a,0) 

Ot? = m , 20 

J 7 (a, 0) 

as a function of temperature both numerically [by the direct use of the Lifshitz formula Q] 
and analytically [using Eq. fTTj) ]. In Fig. 2 the numerical results are presented by the solid 
line and the analytic ones by the dashed line for a dielectric with e D = 1.1 (77 = 0.1) at a 
separation a = 1 fim. As is seen from this figure, at T < 220 K there is practically exact 
coincidence of numerical and analytic results. What this means is that at r < 1.2 (the 
effective temperature here is equal to T e s = 1145 K) one can use the asymptotic expression 
of Eq. (JT7j) in order to get accurate results. 

The asymptotic regime of high temperatures is given by Eq. (JHJ) • This regime is achieved 
at T = 1500 K where the values of StJ 710 , calculated by Eqs. and Q, differ for only 
0.2%. For higher temperatures (r > 8.2) Eq. (J£JJ) is applicable for accurate calculations. 

We now turn to the dielectrics with larger e®. For the sake of convenience in the com- 
parison of experiment and theory, we will present the results of the calculation as a function 
of separation rather than temperature. The point is that in all the previous experiments 
of Refs. the temperature was constant (T = 300 K) and 
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the separation distance was variable (there is only one proposed but not yet completed 
experiment 45, |^ exploiting the effect of two different temperatures). Theoretically, the 
two asymptotic regimes of low (T <C T e g) and high (T 3> T g) temperatures are equivalent 
to the limits of small and large separations, respectively. This becomes evident if it is re- 
called that T e g ~ a~ l . We consider the separation region 100 nm < a < 1400 nm where the 
nonmonotonous behavior of the thermal correction to the Casimir energy is observed for 
dielectrics with larger e (not dilute). In fact, even in that separation region effects which 



are due to t 
Refs 



q ue i ii 



ie nonideality of a metal and the absorption bands in a dielectric (see, e.g., 
]) become important. At shorter separations (a < 100 nm) these effects lead 
to the result that the model of an ideal metal and a dielectric with constant dielectric per- 
mittivity becomes inapplicable (see Sec. IV where the case of real metal and dielectric is 
considered) . 

In Fig. 3a we present the values of the relative thermal correction ()20|) to the Casimir 
energy at T = 300 K as a function of separation computed by the use of Lifshitz formula 
(JHJ) for ideal metal and dielectrics with different dielectric permittivities = 3, 6, 7, 10 
(lines 1, 2, 3, and 4, respectively). As is seen from Fig. 3a, lines 3 and 4 demonstrate the 
nonmonotonous behavior of the thermal correction to the free energy (recall that in the 
case of an ideal metal and dilute dielectric only the Casimir pressure is a nonmonotonous 
function of the temperature, see Fig. 1). Curiously, line 3 possesses both a minimum and a 
maximum values within the separation interval under consideration. 

For line 4 (e$ = 10) the zero value of the thermal correction is achieved at a\ = 0.2 fim 
and at = 1.25 /im. Within the separation region a\ < a < the relative thermal 
correction is negative. Using the respective effective temperatures (T e ^ = 5724 K at a = a x 
and T$ = 916 K at a = 0-2) we find that in terms of the dimensionless parameter r, 
introduced in Sec. II, the thermal correction is negative within the interval 0.33 < r < 2.06. 
The minimum on line 4 is achieved at a = 0.9 /im (r = 1.48). What this means is that 
the entropy of a fluctuating field takes negative values within the interval 0.33 < r < 1.48 
(in the case of a plate made of dilute dielectric the entropy is positive). In terms of usual 
temperatures for the plates at a separation a = 0.2 fim, the Casimir entropy is negative 
within the interval 300 K < T < 1350 K (if, as assumed in this section, e D is independent of 
the temperature). Note that for the dielectric with = 10 the thermal correction is almost 
0.5% at a = 0.9 /mi, i.e., much greater than in the case of a dilute dielectric. It remains, 
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however, still too small to be observed experimentally. 

To illustrate the behavior of the Casimir entropy as a function of temperature in the 
configuration of a plate made of ideal metal and another plate made of dielectric with a 
frequency-independent dielectric permittivity, we plot it in Fig. 3b for = 7 and separation 
distance a = 600 nm. As is seen from Fig. 3b, the entropy is negative within the interval 
from 137K to 311 K (i.e., from 0.45 to 1.02 in terms of r). The minimum value of the entropy 
equal to — 14KeVm _2 K _1 is achieved at a temperature T = 238 K (r = 0.78). When this 
result is compared with the above case of — 10, it is apparent that the region, where the 
entropy is negative, is narrowed with the decrease of the dielectric permittivity. Thus, for 

< 6 the Casimir entropy is already positive at any temperature. 

The negativeness of the Casimir entropy within some temperature (separation) intervals 
should not become of concern. It is self-evident that the entropy of a closed system, which 
includes the space occupied by the dielectric plate, is positive (the second plate is made 
of an ideal metal with the Dirichlet boundary condition on it; for this reason it does not 
contribute to the entropy of the system under consideration). It should be noted also that 
the Nernst heat theorem is satisfied perfectly good. To proof this fact analytically, we apply 
the Abel-Plana formula (jllj) in Eq. Q where the role of F(l) is played by the integral in 
the right-hand side of Eq. (jSJ). Preserving only the lowest expansion order in the parameter 
rl, we arrive at 



F(/)=r 2 / 2 f-j f dy (21) 

fc fc -f 1 



X 



L(s? + l)-( £ ?-l)e- 
From this after integration it follows 



e 



-y 



— + 0(tH 3 ). 

y 



<V D - ir 

F(it) - F(-it) = m¥ l ° L (22) 

Performing the integration with respect to t in Eq. pip, we obtain from Eq. (jSJ) the asymp- 
totic behavior of the Casimir free energy 



^ V T)^V^0)-^(f^f. (23) 

Note that for dilute dielectric this coincides with the term of Eq. (JT7j) of order r 3 as it should 
be. 



2 
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From Eq. (J2~3*j) we finally arrive to the low-temperature asymptotics of the entropy 

which is also in perfect agreement with the second-order term in Eq. ()18j) obtained for dilute 
dielectrics. From Eq. (|24|) it follows that the Casimir entropy goes to zero as the second 
power of the temperature, i.e., the Nernst heat theorem is satisfied. The comparison with 
the numerical computations shows that the asymptotic expressions fl2rij) . (f2~4"|) work good for 
all r < 0.1. 



IV. THERMAL CORRECTIONS TO THE CASIMIR INTERACTION BETWEEN 
REAL METAL AND DIELECTRIC 

In this section we consider one of the plates made of real metal (Au) and the other plate 
made of a real dielectric (Si or a-Al 2 03). Both these chosen dielectrics possess relatively 
large values of the static dielectric permittivity and quite different behavior of £ D (i£) around 
the characteristic frequencies for the separations under consideration. The Casimir free 
energy is found by the use of the complete Lifshitz formula describing the case of two 
parallel plates made of real materials. The dielectric permittivity of real materials along the 
imaginary frequency axis can be obtained through the dispersion relation 

=i + ! r^-^M. (25) 

Here the imaginary part of the dielectric permittivity is calculated as 2ni7i2 where n\ and 
n,2 are the real and imaginary parts of the complex refraction index tabulated, e.g., in 
Ref. 49]. For Au the available tabulated data are extended for lower frequencies using the 



usual procedure (see, e.g. , Refs 



can be found in Refs. 



21 



24 



47 



48]). The resulting behavior of function of £ 

. The same procedure, applied in the case of Si (here 
the tabulated data for lower frequencies, than for Au, are available so that no additional 
extension of data is needed), leads to the results shown by line 1 in Fig. 4. 

There are also good analytic formulas describing the behavior of the dielectric permittiv- 
ities of different materials along the imaginary frequency axis. As an example, the dielectric 
permittivity of a-Al 2 03 is well described j^| in the Ninham-Parsegian representation 

£ «) = 1 + g2 + e , (26) 
13 



where u>m = 1 x 10 14 rad/s, cjuv = 2 x 10 16 rad/s are the characteristic absorption frequencies, 
and Cir = 7.03, C*uv = 2.072 are the corresponding absorption strengths in the infrared 
and ultraviolet ranges, respectively. The dielectric permittivity of a-A^O^ as a function of 
£ is plotted in Fig. 4, line 2. As is seen from Fig. 4, the dielectric permittivities of Si and 
OJ-AI2O3 are qualitatively different in the region of characteristic frequencies £ c ~ 10 15 rad/s. 
In fact, for Cc-Al 2 03 in the region around £ c the values of e D correspond to the second step 
of line 2 and are several times less than = 10.1, whereas for Si the static value of the 
dielectric permittivity = 11.66 is preserved up to the region of characteristic frequencies. 

It should be stressed that the computations below are unaffected by the controvercies 
concerning the contribution of the zero-frequency term of the Lifshitz formula in the case 
of real metals (this contribution is different in the approaches using the Drude dielectric 
function and the Leontovich surface impedance, see Introduction). The reason is that in our 
case only one plate is made of a real metal, whereas the other one is made of a dielectric. 
According to Eq. (J2J, for dielectrics with finite £ D (0) = it follows 

r^(0,y) = £ 4^, r?(0,y) = 0. (27) 

As a result, if one plate is made of a dielectric, the transverse electric mode at zero frequency 
does not contribute to the Casimir free energy regardless of the approach used to describe 
the metal of the other plate [i.e., regardless of the value of the transverse reflection coefficient 

rf(0,y)}. 

A further distinctive feature of the different approaches to the thermal Casimir force in the 
case of real metal, which might play a part in determining the contribution to the free energy 
at nonzero Matsubara frequencies, is the form of reflection coefficients. In the framework of 
the impedance approach fi Q, in p.ace of the nsna, reflection coefficients 0, expressed 
in terms of the dielectric permittivity, one should use the coefficients, expressed in terms of 
the Leontovich surface impedance. This, however, does not present a problem because, as 
was demonstrated in Ref. |l4j, at all nonzero Matsubara frequencies the contributions from 
both types of the reflection coefficients are practically the same. 

In Fig. 5 we present the results of the calculation for the relative thermal correction 
to the Casimir energy between Au and Si plates obtained by Eqs. and (J2J) using the 
procedure described above at T = 300 K (solid line). For comparison, in the same figure the 
dashed line shows the results obtained using the approach of Sec. Ill (i.e., for ideal metal 
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and the dielectric with a frequency- independent dielectric permittivity Eq = 11.66 equal to 
the static permittivity of Si). As is seen from Fig. 5 (solid line), there is a wide separation 
interval 0.2 /mi < a < 1.3 /im where the relative thermal correction to the Casimir energy 
in the case of real materials is negative (in terms of the dimensionless variable this holds for 
0.33 < t < 2.14). The minimum value of the thermal correction 5tT = —0.006 is achieved 
at a = 0.95 /mi (r = 1.56). What this means is that the Casimir entropy in the case of real 
materials is negative within the separation region 0.2 /im < a < 0.95 /mi (or, in terms of r, 
for 0.33 < t < 1.56). The comparison with the dashed line shows that for Si the simple 
model, used in Sec. Ill, leads to the same qualitative results with only minor differences in 
the minimum values of 5tT and the width of the intervals where the thermal correction and 
the Casimir entropy are negative. 

We now turn to the Casimir interaction of Au plate with a plate made of a-Al 2 03. As was 
discussed above, the behavior of the dielectric permittivity of a-Al 2 03 along the imaginary 
frequency axis is different from that of Si. The results of the calculation for the relative 
thermal correction to the Casimir energy as a function of separation at T = 300 K, obtained 
by Eqs. (JTJ), (j2J) and (J2EJ), are shown in Fig. 6 by the solid line. The dashed line is calculated 
for an ideal metal and a dielectric with the frequency-independent dielectric permittivity 

= 10.1 equal to the static dielectric permittivity of a-Al 2 03. As is seen from Fig. 6, in this 
case the solid line presents the monotonously increasing positive function of the separation 
distance. The respective Casimir entropy is also nonnegative within the separation region 
reflected in the figure. The application of the simplified model of Sec. Ill to a-Al 2 03 leads to 
qualitatively different results shown by the dashed line in Fig. 6. This line demonstrates the 
negative thermal correction within the separation region from a\ = 0.25 /im to a 2 = 1.27/tm 
and the negative Casimir entropy within the separations from 0.25 to 0.9 /im. Thus, the use 
of realistic data for the dielectric permittivities of the plates is essential for the final results. 

V. CONCLUSIONS AND DISCUSSION 

In the foregoing we have investigated the thermal corrections to the Casimir interaction 
between metallic and dielectric plates. This was done both analytically (using the ideal- 
ized model of an ideal metal and a dilute dielectric) and numerically (for the ideal metal 
and dielectric with a frequency independent dielectric permittivity, and for real metal and 
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two different dielectrics with dissimilar behavior of their dielectric permittivities along the 
imaginary frequency axis). The main conclusion obtained above is that the pressure and 
the free energy of the Casimir interaction between metal and dielectric plates may be a 
nonmonotonous functions of the temperature within some definite regions. This leads to 
the possibility of negative relative thermal corrections and negative values of entropy of the 
fluctuating field (the latter holds only for dielectrics with sufficiently large dielectric per- 
mittivity). Using the proximity force theorem, one can conclude that the relative thermal 
correction to the Casimir force between a plane metal plate and a spherical dielectric lens 
(the configuration used in many experiments) also can be negative. 

The physical interpretation of the obtained results is based on the fact that both the 
free energy and entropy of the closed system under consideration consist of contributions 
from the plates and from their interaction (in the previous sections the latter were denoted 
as J-~ lD , S lT} or as J 7 , S for real materials). The above conclusions about the possibility of 
a nonmonotonous behavior of the free energy and of the negativeness of the entropy are 
relevant not for the closed system but due to the interaction between its parts. In the 
case of two plates made of an ideal metal with the Dirichlet boundary conditions on their 
surface, there is no penetration of the fluctuating field inside the plates. In this case the 
characteristics of the closed system coincide with those obtained for the interaction between 
the plates. As a result, for ideal metals the free energy of a fluctuating field is a monotonous 
function and the entropy is positive. For two dielectric plates 0, m\ or for one dielectric 
and one metal plate this is not necessarily so. 

It is important to keep in mind that only the interaction parts of the free energy and 
entropy depend on a separation distance: T = !F(a,T), S = S(a,T). This leads to two 
conclusions of considerable significance. The first is that the thermal correction to the 
Casimir force (which is minus the derivative of the free energy with respect to the separation) 
can be negative. The second is that the Nernst heat theorem must be valid separately for the 
contribution to the entropy from the interaction between the plates, so that S(a, 0) = 0, and 
for the entropy of the plates. If this were not the case, i.e., if the equation S(a, 0) = f(a) ^ 
were valid (like in Refs. 25, 2(| 27, 2^ for perfect crystal with no impurities), then the Nernst 
heat theorem for the closed system would be violated as the entropies of the plates do not 
depend on a. Both these conclusions were illustrated above by the example of the Casimir 
interaction between metal and dielectric. 
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To conclude, the nonmonotonous dependence of the Casimir free energy on temperature 
and the negative values of the relative thermal correction (as, for instance, was predicted 



for real metals in the approach of Refs. 



25 



26 



27 



28j and for Si in Ref. [50]) are not in 



themselves excluded thermodynamically. Such behavior for real metals is, however, unlikely 
because there is only a small penetration of the electromagnetic fluctuations at the charac- 



30, 3 



teristic frequencies in the interior of a metal (recall that in the approach of Refs. 
where this property of real metals is taken into account, the free energy is monotonous, and 
the relative thermal correction is positive). The decisive theoretical argument to give pref- 
erence to any approach is, thus, the fulfillment of the Nernst heat theorem for the Casimir 
entropy of a fluctuating field for both perfect crystals and crystals with impurities. In par- 
ticular, if one interacting body is made of a metal and the other of a dielectric, the entropy 
of the fluctuating field vanishes when the temperature goes to zero. 
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FIG. 1: Relative thermal correction to the Casimir pressure as a function of temperature at a 
separation a = 2 fim for two plates, one made of an ideal metal and the other of dilute dielectric 
with = 1.001 (solid line) and = 1.1 (dashed line). 
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FIG. 2: Relative thermal correction to the Casimir energy at a separation a = 1 as a function 
of temperature for two plates, one made of an ideal metal and the other of dilute dielectric with 
= 1.1, computed numerically (solid line) and by the asymptotics of low temperatures (dashed 
line) . 
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FIG. 3: (a) Relative thermal correction to the Casimir energy at T = 300 K as a function of 
separation for two plates, one made of an ideal metal and the other of dielectrics with = 3, 6, 
7, and 10 (lines 1, 2, 3, and 4, respectively); (b) Casimir entropy as a function of temperature for 
the dielectric plate with ej* = 7 at a separation a = 600 nm from a metal plate. 
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FIG. 4: Dielectric permittivity of Si (line 1) and CC-AI2O3 (line 2) along the imaginary frequency 
axis as a function of the logarithm of frequency 
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FIG. 5: Relative thermal correction to the Casimir energy at T = 300 K as a function of separation 
for two plates, one made of real metal (Au) and the other of dielectric (Si) (solid line). The same 
is shown by the dashed line for an ideal metal and dielectric with = 11.67. 
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FIG. 6: Relative thermal correction to the Casimir energy at T = 300 K as a function of separation 
for two plates, one made of real metal (Au) and the other of dielectric (CK-AI2O3) (solid line). The 
same is shown by the dashed line for an ideal metal and dielectric with = 10.1. 
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